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Group field theories whose Feynman diagrams describe 3d gravity with a varying configuration 
of Wilson loop observables and 3d gravity with volume observables at each vertex are defined. The 
volume observables are created by the usual spin network grasping operators which require the 
introduction of vector fields on the group. We then use this to define group field theories that give a 
previously defined spin foam model for fermion fields coupled to gravity, and the simpler "quenched" 
approximation, by using tensor fields on the group. The group field theory naturally includes the 
sum over fermionic loops at each order of the perturbation theory. 

I. INTRODUCTION 

Group field theories (GFTs) provide a method for performing a perturbative sum over 2-complexes with weights 
that give a particular spin foam partition function evaluated on that 2-complex. In this way, GFTs are a proposal 
for solving the problem of triangulation dependence in spin foam models and simultaneously provide a method for 
topology change 0, 0, [|[ . 

An interesting problem in spin foam models for quantum gravity is the inclusion of matter. Current proposals include 
point particles [J] in three dimensions and "strings" || in four dimensions. These are motivated by the description of 
3d particles as conical defects and result in a spin foam model with a deficit angle associated to each edge through 
which the particle propagates. In an appropriate limit, the 3d amplitudes give a non-commutative field theory [6J. 
The formalism can be extended to include particles of arbitrary spin. A group field theory for these particles was 
described in 0] which reproduced Feynman diagrams for both gravity and matter. The effective field theory can also 
be cast in the form of a GFT [8| by perturbing the standard GFT for 3d gravity. 

An alternative approach to adding matter in 3d spin foam models was given in Q . The action for 3d gravity with a 
fermion field was discretised on a triangulation following usual spin foam techniques. This resulted in modified vertex 
amplitudes for the spin foam model that contain an edge in the fundamental representation (representing the fermion) 
which then interacts with the geometry via a grasping operator. Unlike the particle model, the fermion model requires 
a sum over all the vacuum loops of the spin foam - as in lattice gauge theory. In [9( , it was observed that the fermion 
model was not topologically invariant since it now contains local degrees of freedom given by the fermions. A GFT 
would be one way to resolve this dependence on the triangulation. In this paper we construct this GFT and we will 
see that it also naturally includes the sum over fermionic loops that appear in standard quantum field theory. 

To arrive at the fermionic GFT, we first consider Wilson loops and volume operators since these are the building 
blocks of the model. The Wilson loops require an extra argument in the fields to control the propagation of the loops 
and the volume operators will require the use of higher spin fields due to the nature of the operators involved. For 
illustrative purposes, we also include the "quenched" fermion model, which neglects half of the loop configurations, 
as this contains all of the main features of the theory without having to include more complicated configurations of 
loops. We then describe the full GFT for the fermion model and briefly comment on its similarities and differences to 
the GFT for particles coupled to the Ponzano-Regge model. 
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II. WILSON LOOPS 



We start by defining the GFT for 3d Euclidean gravity with Wilson loops. The Feynman diagrams of this field 
theory will give Ponzano-Regge spin foam models with Wilson loop observables 10]. In fact, the GFT will give all 
possible Wilson loops in the spin foam so in itself is unlikely to be physically interesting. However, it will be useful for 
defining the matter coupling. For the Ponzano-Regge model defined on a triangulation A, the Wilson loop observable 
Oyy is defined as the trace in the representation J of the holonomy around a closed loop W in the dual triangulation 
A* 



w 



'EW 



The expectation value of Oyy is then given by 



<0&}pr(A) 



SU{2) 



n dge 



O'U l[6(G f .] 



(1) 



(2) 



The group variables g e * represent SU(2) holonomies along dual edges e* of A, Gy* = rie*e/* 5e* is the ordered product 
of holonomies around a dual face /* of A. Instead of the usual 6j symbol, the spin network vertex amplitudes for the 
model will contain an extra edge in the spin J representation at any vertex that is traversed by the Wilson loop. 




We can now define the GFT for a Wilson loop observable. Since the GFT sums over 2-complexes, we will in fact create 
all possible Wilson loops at each value of the GFT expansion parameter. This will include disjoint loops so from now 
on we use W to denote a configuration of any number of Wilson loops. Unless the loop configuration W completely 
saturates the spin foam, we will also require vertex amplitudes for those vertices that are not part of the loop. For 
this, we recall the pure gravity field 11] defined as a map from three copies of SU(2) to the complex numbers 



0(51,52,53) : SU(2) x SU{2) x 5(7(2) C 
As usual, we require that this is invariant under permutations a of the group variables 

0(51,52,53) = 0(0<r(l), 5<t(2), 9a(3)) 

and we make the field SU (2) invariant 

0(551,552,553) = 0(51,52,53) ; 5 G SU(2) 

with the projector P a 



-Pa0(5l,52,53) 



da 4>(agi,ag 2 ,ag 3 ) 



5(7(2) 



(3) 
(4) 
(5) 

(6) 



With the field defined as above, we can write the action S GR for the GFT whose Feynman diagrams produce Ponzano- 
Regge spin foams. 



'gr[0, A G r 



] = 2 J n d 5i^i0(5i,52,53)-Pa 2 0(5i,52,53) 


Y[ dg t P ai (f>{gi, g 2 , 53)fa 2 0(5i, 55, 56)-Pq 3 0(52, 54, 5e)-Pa 4 0(53, 54, 5s) 



(7) 
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To include the Wilson loops, we define an additional field ip with an extra argument 

^(Si, 32,5355) : SU(2) x 5(7(2) x SU(2) x SU(2) -> C (8) 

This extra variable g will control the way the Wilson loop propagates throughout the 2-complcx. We project the 
additional argument of the field onto the representation J specified by the Wilson loop with the following projector 

^ J {91,92,93] g) = (P J ip)(gi,g2,93;g) = / dh d J x J {gh~ 1 )ip(g ll g 2 ,g3;h) (9) 

JSU(2) 

where dj = 2 J + 1 is the dimension of the irreducible representation J and x' J IS the character. We also demand that 
the field is SU{2) invariant by projecting 

P a ip{.gi,92,g3;9) = / da ip(gia,g 2 a,g 3 a;ga), (10) 

JSU(2) 

this will create the four-valent intertwiners in the Wilson loop vertex amplitudes. Finally, we demand permutation 
symmetry on the first three arguments. 

^(51,52,53:5) = i>(g*(i), 3(7(2), 3 CT (3); g) (n) 

We can use the Peter- Weyl theorem to perform a Fourier decomposition of the fields 

^ J (3l,32,3 3 ;3) = / da / dk ^mikXlkimsksrntkt V d n d n d K d H 

3 „ mi .n,M,a,b J SU(2) J SU{2) 

l<i<4 

x D% ini (gi)D% 2n2 (32)^ 3 „ 3 (5 3 )^ 4 „ 4 (h) 
x Di\ ki (a)Di: 2k2 (a)Di* 3k3 (a)D^ ki (a) 

x D J ab {g)D{ a (h- 1 ) (12) 
Performing the group integrals reduces this to 

P a ip J (gi, 32, 33; 3) = €!*m 3 m 4 V^MA 

l<i<4 

^i ni (9i)^ 2 n 2 (!/2)^3 n 3(53)<„Js)C^ 4 (13) 
Where C^ k ^j S is a four valent SU(2) intertwiner labelled by s and we have redefined the coefficients as 

Tmim2m 3 m4 Vm 1 kim2k2m 3 k 3 m l k l ktk^k^ki \ ± ^ : l 

We define the action as a functional of the two fields to be 

S w [(j), tp J , A GR , A w ] = S GR y>, A G r] + \ [ f[dg l P a tp J (g 1 ,g 2 ,g3;g)P a tp J (gi,g 2 ,g3;g) 

ZJ i=i 

A w f 6 

+ ~r W d 9t dgP ai ip J {gi,g2,g3;g)P a2 ip J (gi, 35, 3e; 3)^3^(32, 54, 36)^4^(33,34, 35) 

i=l 

(15) 

With this action we can now compute the propagators and vertex amplitudes of the theory. 
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A. Feynman amplitudes 



We can now consider the partition function for the theory and express it in terms of the Feynman diagram expansion. 
The expansion will implement a sum over 2-complexes and Wilson loop configurations on each 2-complex. 



r 



sym(r) 



Z[T}. 



(16) 



Here T denotes a particular Feynman diagram in this sum. The number of pure gravity vertices in T is denoted u GR [r], 
the number of Wilson loop vertices is v w [T] and the symmetry factor of the diagram is written as sym(r). Each term 
Z{T] in the sum gives the spin foam partition function evaluated on the 2-complex defined by T. 

The Feynman rules for the GFT are as follows. There are two propagators, the pure gravity propagator and one 
with an extra delta function that determines the path of the Wilson loop. 
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Similarly, there are two vertex amplitudes corresponding to a vertex that is not intersected by a loop and one that is. 
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Note that for a fixed v GB , i> w is the number of dual vertices traversed by the Wilson loops. Special cases include 
v w = which corresponds to a pure partition function without observables and v GR = which contains all possible 
Wilson loops that completely saturate the spin foam. 



III. GEOMETRIC OPERATORS 



The geometric operators we consider here are constructed from polynomials of the discrete frame field used in 
the construction of 3d spin foams. Expectation values of these operators can be computed using the technique of 
generating functionals [12j and result in so called grasping operators acting on the spin network vertex amplitudes. 
In particular, we will consider the volume operator although others can be constructed in a similar fashion. There 
are different ways in which to construct a quantum mechanical volume operator in terms of the frame fields which are 
equivalent classically. These have been studied in detail in [13(. For a tetrahedron, t, we define the volume operator 
Vt to be 

ei,e 2 ,e 3 
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The Xl are the discrete frame fields associated to each edge e of A and carry an SU(2) Lie algebra index. The 
summation is over all triples of edges e\,e-2,e 3 in t that give a non-zero volume. To compute expectation values of 
this type of observable, one introduces su{2) valued sources J e on each edge e and defines the generating functional 
Z[J e ] to be 



Z[Je] = TT / dg e , TT / dX d exp ( i V Tr X e (G e + J e ) ) . 

e,JSU(2) \,Jsu{2) \ ~ J 

The (unnormalised) expectation value of the volume operator can then be computed using 

(V t > OB (A) = TT / dg e * J] / dX e , V euKX^X^Xg exp (i VTrX e G e ] 



(18) 



E 

Vei,e 2 ,e 3 



5J/ 2 5J e ^ 



Z[J e 



Je = 



(19) 



The source derivations insert Lie algebra generators into the appropriate edges which are then contracted with 6ijk- 
Remembering that 6ijk is an intertwiner between three vector representations, then, up to some normalisation factors, 
the resulting spin network vertex amplitudes will be of the form 




The GFT: We begin by observing that the grasping operation cannot be implemented simply as a projector on 
the standard Boulatov GFT as the grasping would be non-local in the fields. We resolve this problem by considering 
vector fields on SU(2). With an appropriate invariance property, this will allow us to construct the graspings by 
contracting with the index of the field. We begin by defining a vector field on three copies of SU(2) 

4>a(9i, 92, g 3 ) ■■ SU(2) x SU(2) x SU(2) - Vi (20) 

Where Vi is the j = 1 SU(2) representation space and a is an index with values 1, 2, 3. In order to create the necessary 
intertwiners, we use the following projector on the field 

P a ,0(f>a(gi, g2, 93) = da d/3 Dl b ((3)(j) b (g 1 a,g2a,g 3 l3a~ 1 f3). (21) 

JSU(2) 

Note that with this projection, the field satisfies the following SU (2) invariance property 

D lb( h )Pa,/3(t)b(gih, g 2 h, g 3 h) = P a ,p<t) a (gi, g 2 , 53)- (22) 

We also require the invariance under permutations a of the group variables by summing over permutations of the field 
arguments 

<f>a(9l, 92, 93) =22^a(ga(l),9<T(2),9<r(3))- (23) 

(T 

The mode expansion of this field is then 

P*,M9l,92,g3) = Cifexm 2 fc 2 m 3 &36 ^ d h ^ d h d U ^m 1% (Sl)£>m 2 n 2 (ff2)-D^ 3 „ 3 (<? 3 ) 

l<i<3 

x / da D^ ki (a)D^ k2 (a)Dll q3 {a- 1 ) f dfi Z>£ {/3)D h {/3)D h (24) 

JSU(2) JSU(2) 



G 



Which reduces to 



Pa,M9i, 92,9a) = £ «fi 3 2m3 y/ d h dh dh d H r C™£ a 

ji ,rrii ,m ,ki,qsa : b 
l<i<3 

xD% ini ( gi )D% 2n2 (g 2 )D% 3n3 (g 3 ), 



where the modes have been redefined as 



Ul02]3 

mim277i3 



^ m\k\ra 2 k 2 m 3 k 3 b kik 2 p 3 P3k 3 b 



(25) 



(26) 



n i tki ? P3 

Ki<3 



With this new field, we define the action for the GFT that has a volume grasping on every vertex amplitude. 
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S'vol^, A VOL ] = 2 f Y\. dgiPa i < t > ( 9l > 92 > 93)Pa 3 <i>(gU 92,9s) 

+ 2 J ndgiPauPiMgi, 32, g 3 )P a2 ,i3 2 <Pb(gi, 32, gz)5ab 



+- 



Y\ dgi Pai,Pi<l>a(9u 92,93)Pa*,pi<k(9l,95,96) P '*3,P3 0c (92 , 34 , 36 ) -Pq 4 <M33 , 34 , 35 ) £ab c (27) 



Feynman rules: The partition function and its Feynman diagram expansion is 

r \"VOL[r] 

Z = / D^ a D^e- SvOL[ *'*°" XvOLl = V V ° L , z\T]. 
J ^ sym(r) 



(28) 



With w V ol[F] giving the number of vertices of T. There are two propagators for the theory, one for the pure gravity 
field and one between the grasped fields. 



9i 

92 ■ 

93 ■ 



9i 

92, 
53 



9l > 

92' 
93' 



>9i 

'92 

• 93 



There is a single vertex that is formed by the interaction of three grasping fields with one pure gravity. The dotted 
lines connect with the vector field propagator while the gravity fields use the usual propagator, note that this means 
the vector field cannot propagate to the gravity field. 



315233 



9i < 
95 ' 

96' 



363234 




Note that due to the form of the vertex, the Feynman diagram expansion will not produce all 2-complexes but 
a subset of them. We ignore this issue as we are only considering SVol as a stepping stone to construction of the 
fcrmionic model. Since we have imposed the permutation symmetry on the field, we will obtain all possible graspings 
on the tetrahedral spin network. These different graspings are described in [l3|. 
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IV. FERMIONS 

Using the techniques above, it will now be possible to write the GFT whose Feynman diagrams give fermion fields 
coupled to 3d spin foam gravity. First, we briefly recall the fermion model and describe how to obtain the massless 
case, for more details refer to Q. 

A. Massless Fermions in 3d gravity 

The fermion model is obtained by discretising the gravity sector as usual and then placing the fermion fields at the 
dual vertices v* of A in the same way as conventional lattice gauge theory. Once the fermionic integration has been 
performed, one is left with a sum over two different sets of loop (±,£2 that saturate all the dual vertices of A. As for 
the Wilson loops, these can consist of more than one loop so we refer to £™ x , C^ 2 as loop configurations containing ni 
loops as £™ ! = {£±, ...£ ni },i — 1,2. For each loop Z, we insert the following into the pure gravity partition function 



A e »D*{g e *)-D*{g e *)A e *-A\ (29) 



Where D^{g) is the representation matrix of g in the fundamental representation and the spinor indices have been 
suppressed. The A e * term describes the area of the triangle / dual to e* through which the fermionic loop I propagates 
and is given by 

A e -= tUK^X^X*. (30) 

e 1 .e 2 Ef 

Here, the summation is over edges e in the triangular face / dual to e* and a 1 denotes the Pauli matrices considered 
as an intertwiner between two spinor and one vector representation space. 

In fact, the spin foam model in Q considered massive fermions and the integration over the Grassmann valued 
fermion fields was performed by using a hopping parameter expansion in the inverse fermion mass to calculate the 
Pfaffian of the Dirac matrix defining the action (equation (43) of Q). The massless case can be obtained directly from 
the Grassmann integration, i.e. by repeated use of equation (46) of 0] without the volume term, and one obtains the 
sum over loop configurations given above. Alternatively, using the hopping parameter expansion, if one considers the 
limit in which the massparameter is small then the highest order term in the expansion dominates. This is the term 
r„ in equation (50) of Q which again gives the sum over loop configurations described above. There are two different 
configurations of loops in the spin foam since it is actually the Pfaffian squared which appears in the path integral. 

Thus the partition function Z F (A) for fermions coupled to 3d gravity can be considered as the expectation value of 
a particular observable for the usual GR partition function 



(31) 



Z F (A) = V TT / dg e , J] f dX e , ( TT FJ TT fA exp( l VTrA e G e ) 

This can be computed with the generating functional techniques as for the volume operators. One can see from pip 
that for each loop £, one obtains a Wilson loop in the spin half representation that is then grasped to the edges of 
each triangle it passes through. Up to symmetry, the different vertex amplitudes for the resulting spin foam model are 
shown in figure[TJ Note that the orientation of the loops is important as the graspings act on the triangle corresponding 
to the point that the fermion enters the tetrahedron rather than the point that it leaves (neglecting the conjugate 
term in Ff). In fact this is a simplification of the full fermion model in which the effect of the second term in (|29l) has 
been neglected but all of the important details are still preserved. Again, in order to properly define the generating 
functional, one should use the wedge variables but the description given above should be clear enough for our purposes. 

1. The quenched model 

In the calculation of a fermionic observable, an approximation can be made that neglects one of the fermionic loops, 
say £2- This would correspond to considering a single symplectic majorana fermion in (9[. There is then only one 



8 




FIG. 1: Spin network diagrams for the different possible vertex amplitudes for the fermion spin foam model. The dashed lines 
are in the spin half representation and the curved lines denote the grasping operators in the spin one representation. 



type of vertex amplitude 




(32) 



We can understand this as describing a fermion (the spin half line) entering a tetrahedron through one triangle, which 
is acted on by a grasping operator, then leaving through a different triangle. For clarity we first describe the GFT 
for this model before moving onto the more complicated configurations possible when there are two fermionic loop 
configurations involved. We will require three different fields to construct the vertex amplitude: 

• The pure gravity field <j> defined in Section Ull and obeying the same symmetries. 

4>(9i, 92,93) ■ SU(2) x SU(2) x SU(2) c 



A Wilson loop field that has been projected to the spin half representation. 

^(91,92, 93; 9) = (PH)(9i,92,93-,9) ■■ SU(2) x SU(2) x SU(2) x SU{2) -> C 
This corresponds to the fermion leaving the tetrahedron. 



(33) 
(34) 



• A final field that contains the grasping operator between the fermion line and the vertex at its source. The field 
has three free spin one indices 



Uo{9u 92,93)9) ■■ SU(2) x 5(7(2) x SU{2) x SU{2) -» Vi ® Vi ®Vi, 

the fourth argument of the field is projected to the spin half representation 

I 1 

Qbc(9i,92,g3;g) = (P 2 €abc)(gi,g2,g3;g) 



(35) 
(36) 
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and the following projection will provide the necessary graspings 

r 3 

JSU(2) 1=1 

x Qef(9ia, g 2 f3ia~ 1 /3 1 , gzfcu' 1 fa; gfaa' 1 ^) (37) 

We also allow the held to be unchanged under permutations a of the hrst three arguments 

£afcc(5l,52,53;5) = £ Qbc (g a{l)l g a{ 2), 9a(3);g) (38) 

This held describes the fermion entering the tetrahedron and grasping the appropriate triangle. 

We can see that the action must now be constructed so that the fields appear in the correct configuration to describe 
a fermion entering and leaving the tetrahedron. The mode expansions for <f>, if)? are given above and the expansion 

for the quenched fermionic field £f bc into modes (Ofr^iWfcsm^J def S ives 

\ r PiP2P3a^abc)\9li 92, 53; 9) — 2^, ?mim2TOjm4 -1/ a ji a j2 a 33 °n2?2a°Ti3g3& L>n 4*ic^m q 2 qmi 

ji ,m.i ,rii,q3,a,b,c,s 
l<i<4 

x^ ini ( 5l )^ 2 „ 2 (52)^3„3(93)^i n4 (5)- (39) 
with the following redefinition of the modes 

phhhh _ \ ^ I f jihhh \ r<j232i r 1 ^ 31 r ,Jlj2j3 ^ t (ac\\ 

ki,P2,P3,P4,d,e,f,t 
l<i<4 

The action for the theory is given by 

3 



S Q [0, V, 6 A Q ] = ]- f f[ d 9* d 9 p a 1 4> l2 (gi,g2,g3- 1 g) P a2 Q bc (gi,g2,g3;g) tabc 
J i=i 

+ \ I n^ F "i^(3i'32,33) P a2 <P(gi,92,g3) 

+ -T- / n^- 9 ^ 3 ( P 0lP2l3 3 ai£,abc)(9l, 92, 93', g) Pa 2 "0 ^ (ffl , 35 , 56 5 5 ) ^Q 3 0(32 , 54 , 56 ) ^4 (flfc 1 04 , 05 ) £«■ 

i=i 

Feynman rules: The partition function for the theory is 

/■ 3 i 1 i \«q[t] 

z= J] ^-^iw^aw^A^^r]. (4i) 



As before, Uq[T] denotes the number of vertices in the diagram. In this case, the Feynman rules are more specific 
about the direction of propagation of the fermions throughout the 2-complex. To represent this, we orient the fermion 
(dashed) lines and demand that the orientation is consistent within each diagram. The propagators are given by 

gl g\ 5i 5i 

92 52 

93 93 

y y 53 .93 

The single vertex and corresponding amplitude of the quenched model is 
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2. The unquenched model 

We now construct the full GFT for massless fermions coupled to 3d gravity. All of the important features of the 
model are contained in the quenched case, all that remains are the additional complications caused by the inclusion of 
two independent loop configurations instead of one. As one can see from the fermionic vertex amplitudes illustrated 
above, we must provide a different interaction vertex for each possible configuration and orientation. We will omit 
some of the detail as it is essentially the same as the quenched case and the complications involved are uncnlightcning. 
We begin by defining the fields 

• The pure gravity field <f> defined in Section HI1 and obeying the same symmetries. 

<Kgi,g2,g 3 ) :SU(2) 3 ^C (42) 

• A Wilson loop field that has been projected to the spin half representation and corresponding to a fermion 
leaving the tetrahedron through a triangle 

^Hgu92,g 3 ;g) = (Ph)(gi, g 2 , gr, g) ■■ SU{2f -> C (43) 

• A field containing two Wilson loops corresponding to the ends of two fermion lines at a triangle 

Mgi,92,gr,gJ) = SU(2) 5 -c (44) 

Both Wilson lines are projected to the spin half representation 

^T^isu9^mg,f)= I dhi dixHgK 1 ) / dhi dix^(fh2 1 )ip2(gi,g2,g3;hi,h 2 ), (45) 

JSU{2) JSU{2) 

and the field has the usual invariance property by projection 

Paip(gi,92,93]9,f) = / da tp(gia,g 2 a,g3a;gaja) (46) 

JSU(2) 

• The fermion field £ a b e (gi, g2, 33! g) described above which relates to a fermion entering a tetrahedron. 

• A tensor field that corresponds to two fermions passing through a triangle, one entering the tetrahedron and one 
leaving 

x ab c(gi,g 2 ,g 3 ;g,f) : SC/(2) 5 ^ Vf 3 , (47) 

1 

The fourth and final arguments are projected to the spin half representation in the same way as i^i ' 2 ■ There 
will be a single grasping operation from the fermion entering the tetrahedron that will be provided by 

J sum t=i 

xxIe}(9i a ^g2(3ia~ 1 f3i,g3(32a~ 1 (3 2 ;g(33a~ 1 P3,fa) (48) 
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• A tensor field that describes two fermions entering the tetrahedron through the same triangle. 

Xabc,def(9i,92,g S ;g,f) ■ SU(2) 5 -> Vf 6 (49) 

Again, the fourth and final arguments are projected to the spin half representation and we create the two grasping 
operators with 

r 3 

(^7.A/)(Si.9S.S3iff./) = / n^^ da ^ P (/3i)^V/ 3 2)^(/33)^ s (7i)^ 1 t (72) J D}j73) 

JSU(2) i=1 

xXp g r, s t«(9i7ia~ 1 7i,52/3ia~ 1 /9i,5372/32" lQ! 72" 1 / 3 2; gfya^ 1 (3 3 , fl^a" 1 ^) 
This projection gives the field the following invariance property 

Dl p {h)Dl q {h)Dl{h)D\ s {h)D\ t {h)D) u ^ 

11 

= {P^ 7i aXaildef)(9U92,93;g,f) (50) 

All of the new fields defined above are also invariant under permutations of the first three indices. We can now write 
the action 



11 i ii ii 



if 3 1 

= 2 n d9idg P "iV^(5l,S2,ff3;5) P a2 £abc(9l, 92, g3l g) ^abc 

1 f 3 

+ 2 J \\ d 9i P ai <i>(gi,92,g3) P a2 (t>{gi,g2,g3) 

, r 3 

+ 2 J [[ d 9i d 9 dh PpidXaic (gi,g2,g3;g,h) Pp>.a 2 Xdef (9i,g2, 93;9,h)€ abc e def 

1 f -r 3 -!- -II II 

+ 2 J [[ d 9idgdh Pi3n J a 1 Xabc,def(9i,g2,g3;g,h) P a2 ipl 2 (g 1 ,g 2 ,g3]g,h)e abc e def 

f TT -11 11 

■ / [[dgidgdh Pfo7 j a 1 Xabldef(gi,92,g3;g,h) P a2 %l) 2 2 ' 2 {g 1 ,g b ,gQ;g,h)P a3 4>{g 2 ,g i ,gQ) P ai 4>(93,9i,gb)^abctdef 
J »=i 

, 6 

-1— 1- 1 1 1 1 

[[dg t Pp iai Xabc (9i,92, 93; 9, h) Pp> ja2 Xdef (9i,95, 5e5 g, h)P a3 (t>{g 2 , 54, go) P ai 4>(93, 9i, gs)tabc£def 

i=l 

A f ^ 

+ Y / Y[d gi dgdh Pp njai xl£ def (gi,g2,g3;9,h) P a2 ip*{gi,g5,96;g)Pa 3 il>*{93,94,95;h) P a4 (j>(g2, 54, ge)e a b c edef 
J i=l 

A f ^ 

+ Y / \\d gi dgdh Pp tai xI b I(9i,g2,g3\g,h)P a2 Q bc {gi,g 5 ,g^g)P a3 ^^ 
J i=l 

A P /■ -A- 11 -1 -1 

+ 1T / 11 dgidgdh P ai ipi ' 2 (51, 52, 535 3, h)P a2 Q bc ( gi ,g 5 , g 6 ; g)P a3 ^ bc (g 2 , 54, 56! g)P a4 <f>(g3, 54, g^abc^def 

•* i=l 

A /* ^ 1 

' J^c«5id5(i/i P ai Q bc (gi,g2,g3\g)P a2 £l bc (gi,g5,g6\h)P ol3 ^{g 2 ,gA,g6\g^^ (51) 



Ap 
' 4 



A F ' 6 



i=i 
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Feynman rules: The partition function for the theory is 

J a,...,m=l r ^ ' 

The four types of propagator are (in the order given in the action) 

h — — — — h h — — 4- — — h 

9 -> 9 9 -> 9 9 > 9 

9i 9i 9i 9i 9i 9i 9i 9i 

92 92 92 92 92 92 92 92 

93 33 93 33 33 33 33 33 

The vertices, listed in the order given in the action, and their corresponding vertex amplitudes are 




The Feynman diagrams for this GFT reproduce the spin foam model for fermions coupled to 3d gravity defined 
in Q with the simplification mentioned previously and the relevant numerical factors absorbed into the coupling 
coefficient. Note that the fermion spin foam model is only defined on a proper triangulation, not any 2-complex. Some 
of the Feynman diagrams of the GFT will therefore produce spin foams that do not relate to the original definition 
of the model. For example, those in which two adjacent tetrahedra are glued on more than one face and the fermions 
propagate though both of these faces. However, in this particular degenerate example, one can show that the amplitude 
is zero. 
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We can also now see that at each value of v F [T] we obtain all possible loop configurations that completely saturate 
the spin foam. Thus the GFT naturally includes the matter diagrams as well as the sum over 2-complexes. This will 
of course include all matter diagrams on all topologies available for a given value of v F [T] . 

V. CONCLUSIONS 

We have altered the standard GFT for 3d gravity to include observables in the form of Wilson loops, volume operators 
and fermions. While the physical significance of summing over Wilson loops and volume operators is unclear, they are 
both essential for constructing the fermionic GFT. 

The fermionic GFT gives a way to deal with the triangulation dependence in the spin foam model for fermions, one 
of the issues mentioned in 0] and, as one would expect, implements the sum over fermionic loops at the same time. 
We now briefly point out the differences between the fermionic GFT and the GFT for point particles. Firstly, the 
corresponding spin foam models have very different vertex amplitudes, even for the spinning particle. The particles 
live on the edges of the spin foam since they have an associated deficit angle but the fermions propagate along the dual 
edges as they make use of the SU(2) variables to parallel transport the spinors. Secondly, the particle model allows 
3 and 4-valent interaction vertices whereas these are not allowed in the fermionic spin foam due to the Grassmann 
integration. In the fermionic model, there are two separate configurations of fermionic loops and the loops must 
necessarily be closed and saturate all of the vertices (again due to the Grassmann variables) whereas there is only a 
single arbitrary graph in the particle model. Modifications could of course be made to the interaction terms in the 
particle model to remove some of these differences and considering fermionic observables would change the allowed 
amplitudes of the fermion model. 

It would be desirable to give an alternative GFT for the fermion model which takes several Grassmann variables as 
its arguments. In this way it may be possible to recreate the model with fewer interaction vertices, however, it would 
also need to include a dependence on the frame fields in a similar way to Note that this suggestion is not the 
same as that in [l5| where the fields themselves are Grassmann valued. 
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